Electronic correlations decimate the ferroelectric polarization of multiferroic HoMn2 05 
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We show that electronic correlations decimate the intrinsic ferroelectric polarization of the recently 
discovered class of multiferroic manganites RMn205, where R is a rare earth element. Such is 
manifest from ab initio bandstructure computations that account for the strong local Coulomb 
interactions between the manganese 3d electrons -the root of magnetism in these materials. When 
including these the computed electronic, magnetic and lattice structure of multiferroic HoMn2 05 
results in an amplitude and direction of polarization that is in accordance with experiment. The 
microscopic mechanism behind the decimation is a near cancellation of the ionic polarization induced 
by ferroelectric lattice displacements and the electronic one caused by valence charge redistributions. 

PACS numbers; 71.45.Gm, 71.10.Ca, 71.10.-w, 73.21.-b 



Introduction Multiferroics, single phase compounds in 
which magnetism and ferroelectrity coexist, are rare 
0]. These materials, such as for instance the manganites 
RMnsOs (R=Ho, Tb, Y, Eu, etc.) 0, i, i, 0, & 1, , 
are currently of great interest because of the possibility 
to control magnetic properties by electric fields and vice 
versa. It was recently discovered that in for instance 
TbMn205 the magneto-electric coupling is so strong, 
that the electric polarization can be reversed by an exter- 
nal magnetic field [1]. This breakthrough can open new 
routes for the design of magneto-electric devices. 

From a fundamental point of view, however, these mul- 
tiferroic manganites contain a puzzle. In regular, non- 
magnetic ferroelectrics the size of the macroscopic polar- 
ization P computed by modern ah initio bandstructure 
methods agrees exceptionally well with the ones observed 
experimentally [ll| . In the multiferroic manganites, how- 
ever, state of the art ah initio computations predict a P of 
around 1200 nCjcTr? (Tb [13,111] and Ho), whereas the 
experimentally observed values are more than an order 
of magnitude smaller (P=45, 65, 100, and 115 nC/cm? 
for Tb, Ho, Y and Eu, respectively 0, H, i, 0, i, i, S ) ■ 
The question arises whether this large discrepancy is due 
experimental artifacts, for instance the formation of fer- 
roelectric domains, or due to an incompleteness in our 
understanding of the physical properties of these mag- 
netic ferroelectric materials. The outcome of this puzzle 
is not only of fundamental interest, as large theoretical 
values of P promise experimentalists a boost of polar- 
ization upon enhanced material quality, increasing the 
multiferroics' application potential. 

We will show in this Letter that the small polarization 
is intrinsic and caused by electronic correlations. It arises 
because the two contributions to P, the ionic part from 
the lattice displacements and the electronic part from the 
valence electrons are opposite and almost canceling each 



other. In this way the electron-electron interactions drive 
a decimation of the resulting net polarization. We com- 
pute P to be in close agreement with the experimental 
value only when the strong local Coulomb interactions 
between the manganese Sd electrons are accounted for. 

Structure of RMn2 O5 In the following we will focus 
on the case R=Ho, but our conclusions are generic for 
this class of compounds. Neutron and X-ray diffraction 
studies show that these manganites have space group 
Pbam but it is expected that in their multiferroic state 
the actual symmetry group is Ph2im, which allows for a 
macroscopic electric polarization along the b axis @, S 0] ■ 
The orthorhombic Pbam crystal structure of HoMn205 
consists of connected Mn^+Oe octahedra and Mn'^+Os 
pyramids (see Fig. [1]). The octahedra share edges and 
form ribbons parallel to the c axis. Adjacent ribbons are 
linked by pairs of corner-sharing pyramids. Below 38 K 
HoMn205 a commensurate magnetic structure develops 
with propagation vector k = (i, Oj_t) and simultaneously 



the system becomes ferroelectric 

We expand upon previous ab initio calculations by in- 
cluding the very strong local Coulomb interactions be- 
tween the manganese 3d electrons -the Hubbard U. We 
use the projector augmented-wave method (PAW) and 
plane wave basis sets as implemented in VASP [l3| ■ Ex- 
change and correlation are treated using the general- 
ized gradient spin-density approximation [l6| (SGGA) 
and the SGGA+U method |l3, [H. We performed the 
SGGA-f-U calculations for U =4 and 8 eV and a Hund's 
rule exchange of Jh = 0.88 eV for the Mn d-electrons, in 
accordance with values obtained from constrained den- 
sity functional calculations on perovskite mangnites [ll| . 

Starting from the experimental centrosymmetric Pbam 
crystal structure we relax unit cell parameters and ionic 
positions both in the SGGA and SGGA+U schemes, al- 
lowing for the lower symmetry Pb2im structure to de- 
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FIG. 1: (Color online) Schematic view of the crystal struc- 
ture of HoMn205 consting of connected Mn*"'"06 octahedra 
(yellow) and Mn^+Os pyramids (red). The magnetic struc- 
ture of the ground state, labelled A, and its enantiomorphic 
counterpart, labelled R, are shown. The white bars connect 
Mn ions of the Mn"^"'"-Mn'*''"-Mn"^"'" structure along the b axis. 



TABLE I: Computed structural parameters of HoMn2 05 in 
Pb2\m crystal structure using SGGA, SGGA-fU. Distances 
are in A; atoms occupying equivalent Wyckoff positions are 
shown only once. 



velop. The atomic positions are relaxed in a 2x1x1 mag- 
netic super cell (containing 64 ions) along the magnetic 
kx direction [2^. We find that the experimental mag- 
netic structure of HoMn2 05 (labeled by A in Fig. [T|) 
with Pb2im symmetry is indeed the magnetic ground- 
state. The calculated structural parameters for U—0 and 
U=8 are shown in Table HI They are in good agreement 
with both the experimental data and first principles elec- 
tronic structure computations on TbMn2 05 without cor- 
relations 0, • The ionic displacements are small 
but significant, in agreement with the fact that experi- 
mentally the low symmetry structure cannot directly be 
determined In Table |TT] we report the computed 
band gap A and energy gain due to the ferroelectric dis- 
tortion SEpE (see Fig.[3|). 

Ferroelectric atomic displacements The the relaxation 
results in Mn^+ and O3 having significant atomic dis- 
placements along the b direction, compared to which the 
displacements for other ions and in other directions are 
small. In Fig. [2] the displacements of these two types of 
atoms are indicated. In the a and c direction the ionic 
displacements are mirror symmetric so that they will not 
contribute to developing a ferroelectric polarization. 

One qualitative difference between the relaxed unit 
cells obtained in SGGA and SGGA-I-U appears in the 
oxygen octahedra surrounding the Mn'*+ ions: SGGA 
calculations show that the Mn^+ move along the b direc- 
tion and become off-centered (the bonds along the longest 
axis of the octahedron become 1.921 and 1.935 A, respec- 
tively) while switching on the Coulomb interaction U re- 
sults in a suppression of this off-centering. The effect is 
not unexpected as the inclusion of the electronic correla- 
tions increases the band gap significantly, rendering the 
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TABLE II: Gap A, energy gain of the ferroelectric state SEfe 
and Born effective charges of different manganese ions in mul- 
tiferroic HoMn205 with Pb2\m symmetry within SGGA and 
SGGA-(-U. 



electronic system more rigid and less susceptible to per- 
turbations. The instability to off-centering is well known 
to be related to the hybridization of the transition metal 
3d states with the oxygen 2p states [H . Increasing U 
leads to a larger gap, a larger splitting between occupied 
oxygen p and empty Mn d states and therefore a smaller 
effective hybridization between the two. 

Along the 6-direction HoMn205 exhibits a charge and 
spin ordering that can schematically be denoted as a 



chain of Mni^+-Mn^^ 



-Mn-;^ 



In the undis- 



, see Fig. H 

torted Pbam structure the distances d^^ (between Mn!^"'' 
and Mn^+) and d^f,- (between MnJ+ and Mn|+) are the 
same. Relaxation reveals a shortening of distances be- 
tween parallel spins Mn?"*" and Mnl"'' ions -in the ferro- 



electric Pb2im structure (iffj_< dij,^, which optimizes the 
double exchange energy [J . 

Ionic and electronic polarization The total polarization 
P for a given material is the sum of the ionic polariza- 
tion Pion and electronic contribution one Peie 



11, m 



In our ab initio calculations the ionic contribution Pion 
is easily obtained by summing the product of the ionic 
displacements from the centrosymmetric to the ferrolec- 
tric structure with the nominal charge of the ions' rigid 
core. To calculate the electronic contribution PeZe we use 
the Berry phase method developed by King-Smith and 
Vanderbilt within the PAW formalism 2^ [23 . 



First we consider the magnetically ordered high sym- 
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FIG. 2: (Color online) Left: arrangement of the ions in the 
unit cell with arrows indicating the ionic displacements of 
Mn^"*" and O3. Right: schematic view of the magnetic and 



charge ordered Mn!^"^-Mn^ 
direction. 
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arrangement along the 6 



metry Pbam phase of HoMn205, in which by definition 
Pion=0. The electronic part to the polarization, how- 
ever, is not required to vanish. In fact, the material is 
bound to have a magnetically induced ferroelectric polar- 
ization as it is a dislocated spin density wave system in 
which the center of symmetry of the magnetic and lattice 
structure do not coincide, providing a basic mechanism 
for multiferroicity l3, 3| • An equivalent point of view is 
that in this oxide symmetry allows for a purely electronic 
part to the magnetostriction [23| . As a consequence 
the spin ordering induces a redistribution of charge on 
crystalographically inequivalent manganese sites. 

We find in the centrosymmetric A structure a result- 
ing polarization of Peie =284 nC/cm^ along the b axis for 
U—0. We checked that inverting all the spins, produc- 
ing the reversed (R) Pbam structure, leads to the same 
polarization in opposite direction. A finite U alters the 
ferroelectric charge redistribution and gives rise to a po- 
larization Peie=-14/-81 nC/cm^ for ;7=4/8 eV in the A 
spin structure 28]. The fact that Peie is induced by the 
magnetic superstructure is immediately clear from a com- 
putation on this system in the ferromagnetic state, in 
which we find all polarization to vanish. 

In the results above it is remarkable that the electronic 
correlation effects induce a sign change of Peie ■ This is a 
real effect caused by changes in electronic structure and 
is not related to geometric constraints in our calculations. 
A polarization flip is possible because symmetry consid- 
erations alone do not fix the sign of the magnetically in- 
duced polarization -the sign of the magneto-electric cou- 
pling. Thus an inversion of the polarization as a function 
of U is symmetry allowed, just as a temperature induced 
sign change of the coupling is possible and indeed ob- 
served in some materials [29|. At the end of this paper 
we present the microscopic mechanism behind this cor- 
relation induced polarization flip. 

In the relaxed Pb2im structure also the ionic con- 
tribution to the polarization comes into play. We 
find P,on=1193/546/576nC/cm2 for C/=0/4/8, concomi- 
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FIG. 3: (Color online) Ionic, electronic and total polarization 
for different values of U in the relaxed Pb2\m structure (left) 
and Peie in the Pbam structure (right). Inset: ferroelectric 
energy gain, SEfe- 



tant with an electronic polarization Pe;e = 12/-287/-493 
nC / cm^ , resulting in a total polarization P=1205/259/82 
nC/cm^ in the magnetic A structure [2^. These values 
are shown in Fig. [3) The value of the polarization for 
[/=8 is in very good agreement with experiment. From 
the computed values of Pion and P^ie we conclude that 
the Hubbard U causes a near cancellation of the elec- 
tronic and ionic contributions to the polarization and ef- 
fectively reduces the polarization in the multiferroic man- 
ganites by over an order of magnitude. 

Origin of the near cancellation The ionic contribution 
to the polarization is driven by fact that (i) the inter- 
atomic Mn distances depend on spin direction (dfi-^)- < 
dfy-^) and {ii) the valence of the two Mn ions approaching 
each other is different, see left panel of Fig. [H Elec- 
tronic correlations reduce Pion by a factor of two, due 
to the increased electronic rigidity which reduces atomic 
displacements, see Fig. [3l In spite of this correlation in- 
duced reduction the computed Pion is still about six times 
larger than the experimental polarization P. 

The electronic polarization Peie arises from a reorga- 
nization of valence charges caused by both ferroelectric 
lattice distortions and magnetic ordering. Both these 
cause changes in covalency, which in turn cause a flow 
of valence electron charge across the material. Therefore 
the effective charge that is displaced by a distortion can 
be much larger than just the bare ionic value. We com- 
puted this Born effective charges of the Mn ions along the 
Mni^'^-Mn^'^-Mn^"^ direction in the ferroelectric Pb2im 
structure, see Table HIl 

When U=0 all Born charges are larger than the nom- 
inal ones, indicative of the distortions inducing appre- 
ciable changes in covalency. The microscopic mechanism 
is that both the Mn^''^ and Mni^"'' ions transfer charge 
to the oxygen atoms connecting them when they move 
closer together, see Fig. |4l In this situation the electrons 
gain kinetic energy because they can hop between Mn!^'*' 
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FIG. 4: (Color online) Schematic view of the two contribu- 
tions to the ferroelectric polarization in HoMn2 05 in the un- 
correlated ((7=0) and strongly correlated limit (large U). In 
the latter the electronic polarization nearly cancels the ionic 
polarization. The labels "Mn''"^^''"'"" indicate Mn ions that 
have a valence of more/less than 3.5+, respectively. 

and Mni^^ without violating the on site Hund's rule [2^ . 
The induced electronic polarizations (Fig. [?]) are opposite 
and cancel each other, in agreement with the very small 
Peie=12 nC/cm^ that we find in Pb2im when U—0. 

When U is large, the situation changes drastically. 
The system becomes more ionic, depleting valence charge 
from the Mn'*+ sites, which approaches a closed shell 
configuration, and increasing it at the Mn^+ sites. Cova- 
lency of the Mn-0 bonds of the former will therefore be 
strongly reduced at the expense of the latter, see Fig. H) 
Indeed we see in Table U that the electronic correlations 
push the Born effective charge of Mn^+ below even its 
nominal value of 4+. 

Conclusions The overall result is that in these strongly 
correlated multiferroic manganites a large electronic po- 
larization develops, which almost as large as the ionic 
polarization, but opposite in direction, see Fig. 3) An 
intrinsically small net polarization of P—82 nC/cm^ re- 
sults, in very good agreement with the experimental 
value. We therefore conclude that electron-electron in- 
teractions decimate the polarization in the multiferroic 
RMn2 05 manganites. Electronic correlation effects are 
thus of prime importance and quantitatively dominate 
the physical properties of these multiferroic transition 
metal compounds. 
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